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Abstract. We study the decay of a prepared state into non-flat continuum. We 
find that the survival probability P(t) might exhibit either stretched-exponential 
or power-law decay, depending on non-universal features of the model. Still there 
is a universal characteristic time to that does not depend on the functional form. 
It is only for a flat continuum that we get a robust exponential decay that is 
insensitive to the nature of the intra-continuum couplings. The analysis highlights 
the co-existence of perturbative and non-perturbative features in the local density 
of states, and the non-linear dependence of 1/tg on the strength of the coupling. 



1. Introduction 

The time relaxation of a quantum-mechanical prepared state into a continuum due to 
some residual interaction is of great interest in many fields of physics. Applications 
can be found in areas as diverse as nuclear [I], atomic and molecular physics [2] 
to quantum information [3] , solid-state physics [4j [5] and quantum chaos [6] . The 
most fundamental measure characterizing the time relaxation process is the so-called 
survival probability P(t), defined as the probability not to decay before time t. 

The study of P(t) goes back to the work of Weisskopf and Wigner [7] regarding 
the decay of a bound state into a flat continuum. They have found that P(t) follows 
an exponential decay P(t) = exp(— t/to), with a rate 1/to which is given by the Fermi 
Golden Rule (FGR), and hence proportional to the effective density of states (DOS) 
for u> = (energy conserving) transitions. 

Following Wigner, many studies have adopted Random Matrix Theory (RMT) 
modeling [8l [9] for the investigation of P(t), highlighting the importance of the 
statistical properties of the spectrum [TO]. Notably in the context of a many- 
particle system, one should understand the role of the whole hierarchy of states 
and associated couplings, ranging from the single-particle levels to the exponentially 
dense spectrum of complicated many-particle excitations [11] , e.g., leading to a decay 
P(t) ~ exp(— Non-uniform couplings also emerge upon quantization of chaotic 
systems where non-universal (semiclassical) features dictate the band-structure of the 
perturbation, leading to a highly non- linear energy spreading [12J . 

Motivation. — Despite all the mounting interest in physical circumstances 
with complex energy landscape, a theoretical investigation of the time relaxation for 
prototypical RMT models is still missing, and also the general (not model specific) 
perspective are lacking. A reasonable starting point for an RMT modeling is the 
characterization of the physical system of interest by a spectral function C{uj) that 
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describes the power spectrum of its fluctuations (the exact definition is given in the 
next section). For an idealized strongly chaotic systems this power spectrum looks 
"flat" , or using an optional terminology taken from different context it is called "white" 
or "Ohmic" . But in more realistic circumstances C(w) is not flat (see some examples 
in [HI [15]), and one wonders what are the consequences. Of particular interest are 
circumstances in which for small frequencies C(lu) oc w 5-1 with s < 1 ("sub-Ohmic" 
spectral function) or s > 1 ("super-Ohmic" spectral function). For such extreme 
non-flatness the conventional Wigner-Weisskopf-FGR picture is not applicable, giving 
zero or infinite rate of decay respectively. For this reason the decay into an s ^ 1 
continuum is the most interesting and challenging case for analysis. 

Scope. — In this paper, we explore a general class of prototype models where the 
initial state decays into a non-flat (sub-Ohmic or super-Ohmic) continuum. We show 
that the survival probability P(t) = g(t/t ) is characterized by a generalized Wigner 
decay time to that depends in a non-linear way on the strength of the coupling. We 
also establish that the scaling function g has distinct universal and non-universal 
features. It is only for the flat continuum of the traditional Wigner model, that we 
get a robust exponential decay that is insensitive to the nature of the intra-continuum 
couplings. In addition to P(t) we investigate other characteristics of the evolving 
wavepacket, namely the variance A-E aprd (i) and the 50% probability width AE caTt ,(t) 
of the energy distribution, that describe universal and non-universal features of its 
decaying component. 



2. Modeling 

We analyze two models whose dynamics is generated by a RMT Hamiltonian 
W = Hq + V, with Wo — diag{i?„} and n e Z. The first one is the Friedrichs model 
(FM) [13], where the distinguished energy level E is coupled to the rest of the levels 
E n =to by a rank two matrix. The second one is the generalized Wigner model (WM) 
[14] , where the perturbation V does not discriminate between the levels, and is given 
by a banded random matrix. In both cases the system is prepared initially in the 
eigenstate corresponding to Eo, and the coupling to the other levels is characterized 
by the spectral function 



C(uj) = -ImlEi 



V Eo+oJ-Ho+M V 



-l 



\V n , Q \ 2 2ir6{u - {E n ~E a )) (1) 

n^O 

where Hq is obtained from Ho by removing the th row and column. An RMT 
averaging over realizations is implicit in the WM case. 

Given a physical system the spectral function C(uj) can be determined numerically 
(see some examples in |12[I15]) and its various features can be understood analytically 
by analyzing the skeleton which is formed by periodic orbits, bouncing orbits and 
taking into account the Lyapunov instability of the motion. In this paper we would 
like to consider the most dramatic possibility of having non-Ohmic spectral function 
which is conventionally modeled as 

C{uj) = 2 7 re 2 |a;| 5 '- 1 e- | " l/aj '= (2) 

The cutoff frequency ui c defines the bandwidth b = gui c of V nm , where g is the density 
of states. In the FM case ±6 is the furthest reachable state (because n^O states are 
not coupled), and therefore the size of the matrix is effectively N = b+1. 




Figure 1. LDoS for the FM and for the WM via direct diagonalization of 
1600 X 1600 matrices with s = 1.5 and e = 1.44. Upper panel: The log-log scale 
emphasizes the universality of the tails up to ui c . Lower panel: The log-linear 
scale emphasizes the difference in the non-universal core component. 



The assumed form Eq.© for the spectral function C(lu) constitutes the natural 
generalization of the standard FM and WM. By integrating Eq.Q over ui we see that 
the perturbation V is bounded provided s > 0. The s = 1 case is what we refer to as 
the flat continuum (Ohmic case), for which it is well known that both models leads 
to the same exponential decay for the survival probability. For s > 2 the effect of 
the continuum can be handled using 1 st order perturbation theory. We focus in the 
< s < 2 regime and consider the s ^ 1 case for which a non-linear version of the 
Wigner decay problem is encountered. 

In the numerical simulations we integrate the Schrodinger equation for c n (t) — 
(n\ip(t)} starting with the initial condition c n = <5„.o at t=0. We use units such 
that q = h = 1, and consider a sharp bandwidth b. The integration is done using a 
self-expanding algorithm [17]. The spreading profile is described by the distribution 
Pt(n) = |c„(t)| 2 , where the averaging is over realizations of the Hamiltonian. The 
survival probability is P(t) — -Pt(O). The energy spreading is characterized by the 
standard deviation AE apld (t) = [^Ji?„- E ) 2 Ptin)] 1 / 2 , by the median E 50 % = Eo, 
and also by the E 2 s% and E 75 % percentiles. The width of the core component is 
defined as AE COIC (t) = E 75% - E 25% . 
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3. Time Scales 



A dimensional analysis predicts the existence of 3 relevant time scales: The Heisenberg 
time t n which is related to the density of states g; the semiclassical (correlation) time 
which is related to the bandwidth uj c ; and the generalized Wigner times to which is 
related to the perturbation strength: 

tn = Zng, t c = 2n/uj c (3) 

*° ~ \T(3-s) sin( S7 r/2) 
where T is the Gamma function. The numerical prefactor that we have incorporated 
into the definition in Eq.Q will be explained later in Section. 7. We shall refer to g^ 1 
and to lj c as the infrared and ultraviolate cutoffs of the theory. Our main interest is 
in the continuum limit. Assuming further that u> c is irrelevant, one expects a decay 
that is determined by the generalized Wigner time to. 

It should be clear that the existence of a cutoff free universal theory in the 
continuum limit for s ^ 1 is not self evident. In fact the natural expectation might 
be to have either infrared or ultraviolate cutoff dependence. Indeed we find that the 
2nd moment of the spreading depends on the w c cutoff, while to is reflected in the 
FM case but not in the WM case. But as far as P(t) is concerned, we find that a 
one-parameter cutoff free universal theory exists. 



\ °' = 1_ 

/ ~ 7o 



4. The LDoS 



Before analyzing the dynamics, it is important to understand the behavior of the Local 
Density of States (LDoS) [14] . which is defined as follows: 

p(u J )=Y / \(v\0)\ 2 ^-(E„-E )) (5) 

where \v) are the eigenstates of the full Hamiltonian H. An RMT averaging over 
realizations is implied in the WM case. Once the LDoS is computed, we can use it to 
calculate the survival probability: 

2 r -i 2 



P(t) 



Mil 



|0) = FT 2tt ( o(w) 



(6) 



where FT denotes the Fourier transform. For flat bandprofile (s — 1), the LDoS 
p(uj) = (l/7o)/(w/7o) is a Lorentzian f(x) — (1/tt)/(1 + x 2 ) [14], leading to a Wigner 
exponential decay for P{t). For (s ^ 1) the ensuing analysis shows that p{oS) has a core- 
tail structure [THl H3 [T^] . Namely, it consists of two distinct regions x ^> 1 and x < 1 
that reflect universal and non-universal features of the problem respectively. The tails 
x 3> 1 can be calculated using 1 st order perturbation theory leading to f(x) cx l/x 3 ~ s . 
This component we regard as universal. The core (x < 1) reflects the non-perturbative 
mixing of the levels, and it is non-universal. In the WM case we argue that for x -C 1 
it is semicircle-like, while for FM we have a singular behaviour f(x) ~ x 1_s . These 
findings are supported by the numerical calculations of Fig.l, and are reflected in the 
behavior of P(t) as confirmed by the numerical simulations of Fig. 2. 



5. Friedrichs model 



Using the Schur complement technique, we can calculate analytically the LDoS for the 
FM. The Green function is G 00 (w) = {[w-A(w)] + i{T(uj)/2)}- 1 with the standard 
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Figure 2. The survival probability P(t) for the FM (top) and for the WM 
(bottom). The time is scaled with respect to to- For all curves in the 
main panels g = 1 and s = 1.5. The WM simulations are presented in log- 
log scale in order to contrast it with the FM results. Inset: further analysis 
displaying Y = —ln[P(t)]/t vs X = t in a log-log plot for representative runs 
with (s,e) = black(0.30,4.43), red(1.00, 3.24), green(1.25, 1.14), blue(1.50, 1.09), 
yellow(1.75, 0.50), showing that the decay in the WM case is described by a 
stretched exponential. The red bold dashed line has zero slope, corresponding 
to simple exponential decay for 8=1. 



notations T{w) = C{uj), 



. lj — lj' 2ir 

= e 2 vr cot {sir/2) M^sgnlw) (7) 

In the last line we performed the limit uj c — > oo (with the limiting expression converging 
in distribution). The LDoS of Eq.© is — (l/7r)Im [Goq(w)] leading to 

_ 1 r M/ 2 
P{UJ) ~ « ( w -AM)2 + (I»/2)2 (8) 
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6. Wigner Model 

The analysis of the LDOS for the WM can be carried out approximately using a 
combination of heuristic and formal methods. Our numerical results reported in Fig. 1 
confirm that the LDoS has l Bt order tails \V n fi/{E n — E$)\ 2 that co-exist with the core 
(non-perturbative) component. We can determine the border 70 between the core and 
the tail simply from the requirement po ~ 1 where 

f°° C(u) doj 

For s > 2 we would have for sufficiently small coupling po <C 1 even if we took the 
limit 70 — > 0. This means that 1 st order perturbation theory is valid as a global 
approximation. But for s < 2 the above equation implies breakdown of 1 st order 
perturbation theory at 70 ~ e 2 /( 2_s ). In the tails T-Lo dominates over V, while in 
the core V dominates. Therefore, as far as the core in concerned, it makes sense 
to diagonalize V with an effective cutoff 70. Following [T5], the result for the LDoS 
lineshape should be semicircle-like, with width given by the expression 

1/2 



AE SC = 



Jo 



(10) 



where above we use the effective bandwidth 70, which replaces the actual bandwidth 
oj c (the latter would be appropriate as in [TS] if we were considering the WM without 
the diagonal energies). The outcome of the integral is AE BC ~ 70, demonstrating 
that our procedure is self- consistent: the core has the same width as implied by the 
breakdown of 1 st order perturbation theory. We note that within this perspective the 
s = 1 Lorentzian is regarded as composed of a semicircle-like core and 1 st order tails. 



7. The survival probability 

In the WM case the function p(w) is smooth with power law tails ~ l/ui 1+a where 
a — 2— s. Thanks to the smoothness the FT does not have power law tails but is 
exponential-like. The similarity with the a-stable Levy distribution suggests that 
P{t) would be similar to a stretched exponential, 

P{t) » exp[-(t/t ) 2 -*} (11) 

The expression for to in Eq.Q is implied by the observation that l/|w| 1+Q tails are 
FT associated with a discontinuity — C|t| Q , where C = [2r(l+a) sin(a7r/2)] _1 . 

In the FM case we observe that the function p{uS) in Eq.® features a crossover 
from w 1 "* for \lo\ <C 70 to T{oj)/uj 2 for \uj\ >• 70. Thus, compared with the WM case, 
the FT has an additional contribution from the singularity at uj—0, and consequently 
by the Tauberian theorem |19j . the survival amplitude has a non-exponential decay, 
that for sufficiently long time is described by a power law: 

2sin((s-l)7r) 



P(t) 



(12) 



(2-s)7T (t/t Q )l- 

The long time behavior is dominated by the non-smooth feature of the core, and not 
by the tails. Comparing the exponential and the power-law we can find the expression 
for the crossover time t' that becomes t' ~ [log \s— l|] 1 ^ 2_,s 'to 3> to close to the Ohmic 
limit (s~l). For s — 1 only the exponential decay survives. We emphasize that the 
cutoff independent behavior appears only after a short transient, i.e. for t > t c . For 




Figure 3. Lower Inset: Examples for the time evolution of A_E,. oro for s=1.5 
and 6=800 in the WM case. Main panel: The extracted departure time versus 
the extracted inverse saturation value. This scatter diagram demonstrates the 
validity of one parameter scaling. Upper Inset: The extracted departure time 
versus the perturbation strength e. The theoretical (dashed) lines are based on the 
to estimate of Eq.Q. The deviations of the departure time from the theoretical 
expectation diminish in the limit ui c — ¥ oo. The o corresponds to b = 400, the * 
to b = 800, and the o to b = 1600. 



completeness we note that for the FM with s=2 we get P(t) ~ | log(t / t' c )\ 2 , that holds 
for t c <t <t' c where t' c = t^ 1 /^ \ For s>2 there is an immediate but only partial 
decay that saturates at the value P(t) — |1— po\ 2 for t > t c . 

8. Spreading 

The distinction between core and tail components becomes physically transparent 
once we analyze the time dependent energy spreading of the wavepacket. Using the 
same time dependent analysis as in the s = 1 case of Ref. [T7] , it is straightforward to 
show that the rise of AE COIQ (t) is at t ~ to, and its saturation value is ~ 7o- Thus 
Ai? corc should exhibit one parameter scaling with respect to to. In Fig. 3 we present the 
results of the numerical analysis. Our data, indicate that the expected one-parameter 
scaling is obeyed. We have verified that the slight deviation (shown in the inset) from 
the expected e dependence is an artifact due to having finite (rather then infinite) 
bandwidth in the numerical simulation. 

The physics of Ai? sprd is quite different and not necessarily universal, because the 
second moment is dominated by the tails, and hence likely to depend on the cutoff 
uj c and diverge in the limit uj c — \ oo. Indeed in the WM case we can use the Linear 
Response Theory (LRT) result of pJl [12] 

AE Bpid (t) = [2(C(0) - C(i))] V2 (13) 

where C(t) is the inverse FT of C{w). This gives the saturated value (2w*e 2 /s) 1 / 2 
as soon as t > t c . We now turn to the FM case. The solution of the Schrodinger 
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Figure 4. Scaled spread AE spr d/(L^^e 2 /s) 1 ^ 2 versus scaled time ui c t for the FM 
and the WM. The linear response theory (LRT) prediction Eq.JTsJ, for the WM, 
as well as the exact result Eq, lll4H for the FM are plotted for comparison. 



equation for c n (t) is well known [2], and (setting Eq=0) can be expressed using the 
real amplitude c(t)=co(t). In particular P{t) = |c(t)| 2 and also the energy spreading 
can be computed in a closed form, with the end result 

A£ Bprd (i) = [(l+c 2 (i))C(0) - c{tf + 2c(t)c(i)] V2 

For t < to we can use the estimates c(t) ~ 1 and c(t) « and c(t) ss 
that A^ Bprd (i) behaves as in Eq. ((T3|) . But for t > to we get 

1/2 



(14) 

-C(t) to conclude 



AE SVId (t) a (l+P(t))C(0) 



(15) 



leading to a saturation value smaller by factor reflecting the non-stationary decay 
of the fluctuations as a function of time. More interestingly Eq. ([T4"|) contains a cutoff 
independent term that reflects the universal time scale to- The numerical results in 
Fig. 4 confirm the validity of the above expressions. We note that in the FM case the 
effect of recurrences is more pronounced, because they are better synchronized: all 
the out-in-out traffic goes exclusively through the initial state. 



9. Summary and Discussion 



In this work we have compared two models that have the same spectral properties, 
but still different underlying dynamics. One of them has an integrable dynamics (FM) 
while the other is an RMT type (WM). This is complementary to our past work [20] 
where we have contrasted a physical model with its RMT counterpart. 

Non-Ohmic coupling to the continuum emerges in various frameworks in physics. 
The general WM analysis might be motivated by the study of quantized chaotic 
systems that exhibit non-Ohmic fluctuations due to semi-classically implied long time 
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power-law correlations. In fact typical power spectra are in general not like "white 
noise" (e.g. [THl [201 H2] ) ■ The general FM analysis might be motivated by studies of 
bound states that are embedded in the continuum as in the single-level Fano- Anderson 
model, with diverse realizations in the molecular / atomic / electronic context and also 
with implication regarding photonic lattices: see [21] and further references therein. 

It should be clear that by considering two special models, we do not cover the full 
range of possibilities: In realistic circumstances the perturbation might have any rank, 
and there might be non-trivial correlations between off-diagonal elements (which was 
in fact the case in 20 J. Still our results, since they relate to two extreme limiting 
models (FM,WM), serve to illuminate the limitations on the universality of Wigner's 
theory. 

In the non-Ohmic decay problem that we have considered a universal generalized 
Wigner time scale has emerged. It is not this time scale but rather the functional 
form of the decay that reflects the non-universality. We find that for "non-Ohmic 
chaos" (WM case) the survival probability becomes a stretched exponential beyond 
the Wigner time scale, which is both surprising and interesting. This is contrasted 
with the "integrable" power-law decay that takes over in the long time limit (FM case) , 
and obviously very different from the Ohmic exponential result. Only the standard 
case of flat (Ohmic) bandprofile is fully universal. 

It is worth mentioning that in a bosonic second quantized language the decay 
of the probability can be re-interpreted as the decay of the site occupation n. If the 
interaction between the bosons is neglected this reduction is exact and merely requires 
an appropriate dictionary. In the latter context each level becomes a bosonic site 
which is formally like an harmonic oscillator, and hence the initially empty continuum 
is regarded as a zero temperature bath. Consequently the decay problem is formally 
re-interpreted as a quantum dissipation problem with an Ohmic (s=l) or non-Ohmic 
(s^l) bath. The time scale to is associated with the damped motion of the generalized 
coordinate n. Optionally P(t) could be related to dephasing, and in this case to is 
reinterpreted as the coherence time, as in Landau's Fermi liquid theory. 
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